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The number of atoms trapped within the mode of an optical cavity is determined in real time 
by monitoring the transmission of a weak probe beam. Continuous observation of atom number 
is accomplished in the strong coupling regime of cavity quantum electrodynamics and functions in 
concert with a cooling scheme for radial atomic motion. The probe transmission exhibits sudden 
steps from one plateau to the next in response to the time evolution of the intracavity atom number, 
from iV>3to7V = 2^1^0, with some trapping events lasting over 1 second. 



Cavity quantum electrodynamics (QED) provides a 
setting in which atoms interact predominantly with light 
in a single mode of an electromagnetic resonator of high 
quality factor Q 0. Not only can the light from this 
mode be collected with high efficiency but as well 
the associated rate of optical information for determin- 
ing atomic position can greatly exceed the rate of free- 
space fluorescent decay employed for conventional imag- 
ing 3\. Moreover the regime of strong coupling, in which 
coherent quantum interactions between atoms and cavity 
field dominate dissipation, offers a unique setting for the 
study of open quantum systems . Dynamical processes 
enabled by strong coupling in cavity QED provide pow- 
erful tools in the emerging field of quantum information 
science (QIS), including for the implementation of quan- 
tum computation [5jand for the realization of distributed 
quantum networks |fj,LD|. 

With these prospects in mind, experiments in cavity 
QED have made made great strides in trapping single 
atoms in the regime of strong coupling jjj, la, IsL llfij - 
However, many protocols in QIS require multiple atoms 
to be trapped within the same cavity, with "quantum 
wiring" between internal states of the various atoms ac- 
complished by way of strong coupling to the cavity field 
\li HH LU4 lii! • Clearly the experimental ability to deter- 
mine the number of trapped atoms coupled to a cavity is 
a critical first step toward the realization of diverse goals 
in QIS. Experimental efforts to combine ion trap tech- 
nology with cavity QED are promising , but have not 
yet reached the regime of strong coupling. 

In this Letter we report measurements in which the 
number of atoms trapped inside an optical cavity is ob- 
served in real time. After initial loading of the intracav- 
ity dipole trap with iV « 5 atoms, the decay of atom 
number A^>3— > 2 ^ 1 ^ is monitored by way 
of changes in the transmission of a near-resonant probe 
beam, with the transmitted light exhibiting a cascade of 
"stairsteps" as successive atoms leave the trap. After the 
probabilistic loading stage, the time required for the de- 
termination of a particular atom number N = 1,2,3 is 
much shorter than the mean interval over which the N 
atoms are trapped. Hence, a precise number of intracav- 




FIG. 1: (color online) Schematic of our experiment. Cs atoms 
are loaded into an intracavity FORT (£f) by way of the trans- 
verse cooling field Q3 and the cavity probe field £4. The trans- 
mitted £4 field is directed to a heterodyne detector (HD), 
allowing real-time determination of intracavity atom number. 



ity atoms can be prepared for experiments in QIS, for 
which the timescales (g _1 w 10~ 8 s) « (t « 3 s), where 
t is the atomic trapping time [i| and Hg is the atom- 
field interaction energy. In the present case, the atom 
number is restricted to N < 3, but the novel detection 
scheme that we have developed may enable extensions to 
moderately larger atom numbers N < 10. 

As illustrated in Fig. ^ our experiment combines 
laser cooling, state-insensitive trapping, and strong cou- 
pling in cavity QED, as were initially achieved in Ref. 

A cloud of Cs atoms is released from a magneto- 
optical trap (MOT) several mm above the cavity, which 
is formed by the reflective surfaces of mirrors (Mi, M2) ■ 
Several atoms are cooled and loaded into an intracavity 
far-off-resonance trap (FORT) and are thereby strongly 
coupled to a single mode of the cavity. The maximum 
single-photon Rabi frequency 2g for one atom is given 
by 5 , o/2ti' = 24 MHz, and is based on the reduced dipole 
moment for the 6S 1/2 ,F = 4 -> 6P 3/2 ,F' = 4' tran- 
sition of the D2 line in Cs at Ao = 852.4 nm. Decay 
rates for the QP3/2 atomic excited states and the cav- 
ity mode at ujq = 2irc/Ao are r y/2n — 2.6 MHz and 
k/2tt = 4.2 MHz, respectively. The fact that go ;§> (k, 7) 
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places our system in the strong coupling regime of cav- 
ity QED 0, giving critical atom and photon numbers 
n = 77(2560 « 0.0057, N Q = 2« 7 / 5o 2 » 0.037. 

The cavity is independently stabilized and tuned such 
that it supports TEMqq modes simultaneously resonant 
with both the F = 4 — * F' = 4' atomic transition at 
Ao and our FORT laser at A^ = 935.6 nm, giving a 
length Iq = 42.2 fim. A weak probe laser £4 excites 
the cavity mode at Ao with the cavity output directed 
to detector HD, while a much stronger trapping laser £p 
drives the mode at \f- In addition, the region between 
the cavity mirrors is illuminated by two orthogonal pairs 
of counter-propagating cooling beams in the transverse 
plane (denoted O3). Atoms arriving in the region of the 
cavity mode are exposed to the (£4, £p, ^3) fields contin- 
uously, with a fraction of the atoms cooled and loaded 
into the FORT by the combined actions of the £4 and 
3 fields 0. For all measurements, the cavity detuning 
from the 4 — ► 4' atomic resonance is Ac = 0. The de- 
tuning of the £4 probe from the atom-cavity resonance is 
A4 = +4 MHz, and its intensity is set such that the mean 
intracavity photon number n — 0.02 with no atoms in the 
cavity. The detuning of the transverse cooling field is 
A 3 = +25 MHz from the F = 3 -> F' = 3' resonance, 
and its intensity is about I3 w 4 x 10 1 mW/cm 2 . 

The field £f that drives the standing-wave, intra- 
cavity FORT is linearly polarized, resulting in nearly 
equal ac-Stark shifts for all Zeeman sublevels of the 
F = 3,4 hyperfine ground states of the §Si/ 2 manifold 
[l5j | . The peak value of the trapping potential is — Uo/h = 
-47 MHz, giving a trap depth U /k B = 2.2 mK . A crit- 
ical characteristic of the FORT is that all states within 
the 6P3/2 excited manifold likewise experience a trap ping 
shift of roughly -U (to within w ±15%) HHEllli, 
which enables continuous monitoring of trapped atoms 
in our cavity and avoids certain heating effects. 

Figure E[a) displays a typical record of the heterodyne 
current i(t) resulting from one instance of FORT loading. 
Here, the current i(t) is referenced to the amplitude of 
the intracavity field | (a) | by way of the known propaga- 
tion and detection efficiencies. The initial sharp drop in 
\{a(t))\ around t = results from atoms that are cooled 
and loaded into the FORT by the combined action of 
the (£4,^3) fields 0. Falling atoms are not exposed to 
£4 until they reach the cavity mode, presumably lead- 
ing to efficient trap loading for atoms that arrive at a 
region of overlap between the standing waves at (Ao, Xf) 
for the (£4, £f) fields. Trap loading always occurs within 
a ±10 ms window around t — 0.025 s (relative to t = 
in Fig. Ufa)). This interval is determined using separate 
measurements of the arrival time distribution of freely 
falling atoms in the absence of the FORT 0, 0] . 

Subsequent to this loading phase, a number of remark- 
able features are apparent in the trace of Fig. GJa) , and 
are consistently present in all the data. The most no- 
table characteristic is the fact that the transmission ver- 
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FIG. 2: (color online) (a) A typical detection record in which 
several (N > 4) atoms are loaded into the trap. Heterodyne 
detection bandwidth is 1 kHz. (b) Histogram of 500 such 
traces, binned with respect to the heterodyne signal \{a)\. A 
digital low-pass filter of bandwidth 100 Hz is applied to each 
trace prior to the computation of the histogram, (c) Com- 
parison of the model prediction for p$ (y = 0.5) (+) with 
the measured positions of the histogram peaks in (b) (O) 
Also shown (x 
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sus time consists of a series of flat "plateaus" in which the 
field amplitude is stable on long timescales. Additionally, 
these plateaus reappear at nearly the same heights in all 
repeated trials of the experiment, as is clearly evidenced 
by the histogram of Fig. |2 (b). We hypothesize that 
each of these plateaus represents a different number N of 
trapped atoms coupled to the cavity mode, as indicated 
by the arrows in Fig. [21 

Consider first the one-atom case, which unexpect- 
edly exhibits relatively large transmission and small 
variance. For fixed drive £4, the intracavity field 
is a function of the coupling parameter g^'^'( r) = 
9oGij sin(feo^) exp(— 2p 2 /w^) where p is the transverse 
distance from the cavity axis (z), fco = 2tt/\q, and Gij 
relates to the Clebsch-Gordan coefficient for particular 
initial and final states (i, f) within the F — 4, F' = 4' 
manifolds [20I El) . Variations in g as a function of the 
atom's position r and internal state might reasonably be 
expected to lead to large variations of the intracavity 
field, both as a function of time and from atom to atom. 

However, one atom in the cavity produces a reason- 
ably well-defined intracavity intensity / oc \(a)\ 2 due 
to the interplay of two effects. The first is that for 
small probe detunings A4, the intracavity intensity Ji 
for one atom is suppressed by a factor / relative to 
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the empty-cavity intensity Iq, where for weak excitation, 
/ w 4Cf » 1 with C\ = q 2 /2kj. A persistent, strongly 
reduced transmission thereby results, since the condition 
[C^'i^(r)] 2 ^> x is robust to large fluctuations in atomic 
position r and internal state. The second effect is that 
the F = 4 <-> F' = 4' transition cannot be approximated 
by a closed, two-level system, since the F' = 4' excited 
states decay to both F — 3,4 hyperfine ground levels. As 
illustrated in Fig. ^b), an atom thus spends a fraction q 
of its time in the cavity QED manifold (4, 4'), and a frac- 
tion p w 1 — q in the (3, 3') manifold. In this latter case, 
the effective coupling is negligible (C^ « 4 x 10~ 4 ), 
leading to an intensity that approximates Iq. Hence, the 
intracavity intensity as a function of time I(t) should 
have the character of a random telegraph signal switch- 
ing between levels (Io,h), with dwell times determined 
by (£ 4 , fl 3 ), which in turn set p j^- Since (£4, O3) drive 
their respective transitions near saturation, the timescale 
Tp ~ 1 /is for optical pumping from one manifold to an- 
other is much faster than the inverse detection bandwidth 
(l/2irB) « 160 /is. The fast modulation of the intracav- 
ity intensity due to optical pumping processes thereby 
gives rise to an average detected signal corresponding to 
intensity I\ « plo + ql\ ~ pig for I\ «C Iq ■ 

This explanation for the case of 1 atom can be extended 
to N intracavity atoms to provide a simple model for the 
"stairsteps" evidenced in Fig. Ufa). For N atoms, the 
intracavity intensity should again take the form of a ran- 
dom telegraph signal, now switching between the levels 
(Jo, Ik), with high transmission Iq during intervals when 
all N atoms happen to be pumped into the (3, 3') man- 
ifold, and with low transmission Ik < I\ anytime that 
1 < k < N atoms reside in the (4, 4') manifold, where 
Ik ~ Ii/k 2 for weak excitation with Ac = A4 = 0. 
The intracavity intensities {Ik} determine the transition 
rates {jk^k-i} between states with k and k — 1 atoms 
in the (4,4') manifold, while determines {jk-i^k} 
for k — 1 — > k via transitions from the (3, 3') mani- 
fold. For the hierarchy of states k = 0,1, ... ,N with 
transition rates {7fc-^fc-i.7fc-i^fc}, it is straightforward 
to determine the steady-state populations p k ■ With 
the physically motivated assignments 7fc-i^fe = 70^1 
independent of k and "fk^k-i = li^o/k 2 correspond- 
ing to I k ~ h/k 2 , we find that p[ N) = 1/ £f =0 (fe!) V, 
where y = 70^1/71^0- Hence, for <C Iq, the pre- 
diction for the average intensity is In ~ Pq N ^Iq, which 
leads to a sequence of plateaus of increasing heights 
In+i —* In In-i as successive atoms are lost from 
the trap 

Figure |2fc) compares the prediction of this simple 
model with the measured values of peak positions in (b). 
The only adjustable parameter is the value y — 0.5, re- 
sulting in reasonable correspondence between the model 
and the measurements. Also shown are values p N ^ for 
y = 0.1 to indicate that it might be possible to enhance 
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FIG. 3: (color online) Histogram of 500 traces such as the one 
in Fig. Ufa) 1 binned with respect to both signal strength | (a) | 
and time t. The signals are filtered first as in Fig. |2fb). 



the resolution for a particular range of atom number by 
framing a given few values N\ , N\ ± 1 in the transition 
region Pq Ni ^ ~ 0.5, where Ni ~ 6 in (c). This could be 
accomplished by adjusting the relative strengths of the 
(£4, fields and hence y. 

Although our simple model accounts for the qualitative 
features in Fig. a quantitative description requires a 
considerably more complex analysis based upon the full 
master equation for N intracavity atoms, including the 
multiple Zeeman states and atomic motion through the 
polarization gradients of the Ife beams. We have made 
initial efforts in this direction |23J for one atom, and are 
working to extend the treatment to N > 2 atoms. 

Beyond these considerations, additional evidence that 
the plateaus in Fig. |3 correspond to definite atom num- 
bers is provided by Fig. [3] Here, the data recorded 
for the probe transmission have been binned not only 
with respect to the value of \(a)\ as in Fig. Efb), but 
also as a function of time. Definite plateaus for \(a)\ are 
again apparent, but now their characteristic time evo- 
lution can be determined. The critical feature of this 
plot is that the plateaus lying at higher values of \(a)\ 
correspond to times later in the trapping interval, in 
agreement with the expectation that N should always 
decrease with time beyond the small window of trap load- 
ing around t = 0.025 s. This average characteristic of the 
entire data set supports our hypothesis that the plateaus 
in I (a) I correspond to definite intracavity atom numbers 
N, as indicated in Figs. [3and|21 Moreover, none of the 
500 traces in the data set includes a downward step in 
transmission after the initial trap loading. 

To examine the dynamics of the trap loss more quan- 
titatively, we consider each atom number individually 
by integrating the "plateau" regions along the |(a(t))| 
axis for each time t. The dashed horizontal lines in 
Fig. 13 indicate the boundaries chosen to define the lim- 
its of integration for each value of N. We thereby obtain 
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FIG. 4: (color online) (a) Experimental results for the time 
evolution of the iV-atom populations $jv(t), which are nor- 
malized such that their sum is approximately unity through- 
out the interval shown, (b) The results of a simple model 
calculation Pjv(£) are fit to the data 4>jv(t) with one free pa- 
rameter, the single atom decay rate F. 



time-dependent "populations" $Ar(t) for N = 0, 1, 2, and 
$>3(*) = Ei^^W. which are plotted in Fig. [fla). 
To isolate the decay dynamics from those of trap loading, 
we plot the data beginning at to = 0.034 s with respect 
to the origin in Figs. EJa) and|3] The qualitative behav- 
ior of these populations is sensible, since almost all trials 
begin with N >3, eventually decaying to N = 2, 1, 0. 

The quantities $Ar(i) are approximately proportional 
to the fraction of experimental trials in which ./V atoms 
were trapped at time t, so long as the characteristic dura- 
tion Atjsr of each plateau far exceeds the time resolution 
of the detection. If the bandwidth is too low, transient 
steps no longer represent a negligible fraction of the data, 
as is the case for transitions between the shortest-lived 
levels (e.g., N = 3 — > 2). We estimate that this ambigu- 
ity causes uncertainties in <J>jv at the 5 — 10% level. 

Also shown in panel (b) of Fig. 0] is the result of a 
simple birth-death model for predicting the time evolu- 
tion of the populations, namely Pjv(i) = — r(iVP/v(t) — 
(N + l)P/v+i(£)), where Pj\r(t) represents the probabil- 
ity of N atoms in the trap. The main assumption of 
the model is that there is one characteristic decay rate 
r for trapped atoms, and that each atom leaves the 
trap independently of all others. Initial conditions for 
N = 0, 1 and 2 for the solution presented in Fig. Efb) 
are obtained directly from the experimental data after 
trap loading, $jv(*o)- Since the plateaus for higher val- 
ues of N are not well resolved, we use a Poisson distri- 
bution for N > 3. The mean /i = 5.2 is obtained by 
solving X!iv=3 e ^^fJ- N /Nl = <&N>3(to)- Given these ini- 
tial conditions, we perform a least-squares fit of the set 



of analytic solutions {P/v(i)} to the set of experimental 
curves {$jv(£)} with T the only free parameter, resulting 
in the curves in Fig. 0Jb) with T — 8.5 s" 1 . Although 
there is reasonable correspondence between Figs. 01 (a) 
and (b), 4>jv(t) evolves more rapidly than does P/v(t) at 
early times, and yet the data decay more slowly at long 
times. This suggests that there might be more than one 
timescale involved, possibly due to an inhomogeneity of 
decay rates from atom to atom or to a dependence of the 
decay rate on N. We have observed non-exponential de- 
cay behavior in other measurements of single-atom trap 
lifetimes, and are working to understand the underlying 
trap dynamics. 

Our experiment represents a new method for the real- 
time determination of the number of atoms trapped and 
strongly coupled to an optical cavity. We emphasize that 
an exact number N = 1 to 3 coupled atoms can be pre- 
pared in our cavity within rj 200 ms from the release of 
the MOT. Although the trap loading is not deterministic, 
N can be measured quickly compared to the subsequent 
trapping time t«3s||. These new capabilities are im- 
portant for the realization of various protocols in quan- 
tum information science, including probabilistic proto- 
cols for entangling multiple atoms in a cavity [ill fl2L fl3| . 
Although our current investigation has centered on the 
case of small N < 3, there are reasonable prospects to 
extend our technique to higher values N < 10 as, for 
example, by way of the strategy illustrated in Fig. Efc). 
Moreover, the rate at which we acquire information about 
N can be substantially increased from the current value 
re|(a)| 2 ~ 10 5 /s toward the maximum rate for optical in- 
formation g 2 1 k > 10 8 /s, which can be much greater than 
the rate for fluorescent imaging set by 7 0- 
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